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©Temperature-Corrected Model of Turbulence in Hot Jet Flows 

A standard turbulence model is corrected for total-temperature gradient and compressibility. 

Langley Research Center, Hampton, Virginia 


An improved correction has been de- 
veloped to increase the accuracy with 
which certain formulations of computa- 
tional fluid dynamics predict mixing in 
shear layers of hot jet flows. The CFD 
formulations in question are those de- 
rived from the Reynolds-averaged 
Navier-Stokes equations closed by means 
of a two-equation model of turbulence, 
known as the /e — e model, wherein effects 
of turbulence are summarized by means 
of an eddy viscosity. The need for a cor- 
rection arises because it is well known 
among specialists in CFD that two-equa- 
tion turbulence models, which were de- 
veloped and calibrated for room-tem- 
perature, low Mach-number, 
plane-mixing-layer flows, underpredict 
mixing in shear layers of hot jet flows. 
The present correction represents an at- 
tempt to account for increased mixing 
that takes place in jet flows characterized 


by high gradients of total temperature. 
This correction also incorporates a com- 
monly accepted, previously developed 
correction for the effect of compressibil- 
ity on mixing. 

One of the two equations of the k-e 
model is 

Ft = P 

where p t is the eddy viscosity, p is the 
mass density, k is the time-averaged ki- 
netic-energy density associated with 
the local fluctuating (turbulent) com- 
ponent of flow, £ is the time-averaged 
rate of dissipation of the turbulent-ki- 
netic-energy density, and Cjj is the sub- 
ject of the present correction, as de- 
scribed next. 

In the uncorrected k—£ model, Cj, has 
the constant value of 0.09. The present 
correction alters the value of C v to ap- 
proximate the effects of the temperature 



The Center-Line Mach Number of a supersonic jet from an axisymmetric nozzle with a plenum total 
temperature of 2,009 R (=1,116 K) and an exit diameter (D) of 3.60 in. (9.144 cm) was measured and 
calculated as a function of axial position (x). 


gradient and compressibility on the 
eddy viscosity. Before presenting the cor- 
rection, it is necessary to define some al- 
gebraic terms: 

• The temperature correction enters 
through a function of the gradient of 
the total temperature normalized by 
the local turbulence length scale. This 
function is given by 

Tg =|VT t |A 3/2 /eT t , 

where T t is the total temperature. 

• The turbulence Mach number is given 
by 

M r =(2 k) 1/2 /a, 

where a is the local speed of sound. 

• The compressibility correction enters 
through a function of the turbulence 
Mach number. This function is given 
by 

f(M t ) = (M? -Ml 0 )H(M r -M t0 ), 

where H(x) is the Heaviside function 
of x (the unit step function of x, 
which is 0 for negative x and 1 for pos- 
itive x), and M r0 is a threshold M z 
value (initially set at 0.1) below which 
it is deemed unnecessary to apply the 
compressibility correction. 

Then the corrected value of C v is 

rj-\S 

1 + 

0.041 + /(M t ) 

It should be noted that in the case of 
zero total-temperature gradient, the cor- 
rected value of Cj, reverts to the prior 
constant value of 0.09. 

The present correction was tested on 
experimental data, in comparison with 
four prior standard corrections to the k- 
£ model. The figure presents an exam- 
ple showing that predictions by use of 
the present correction were in closest 
agreement with the experimental data. 

This work was done by Khaled S. Abdol- 
Hamid and S. Paul Pao of Langley Research 
Center; Steven J. Massey of Eagle Aeronau- 
tics, Inc.; and Alaa Elmiligui of Analytical 
Services & Materials, Inc. Further informa- 
tion is contained in a TSP (see page 1 ). 
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C„ = 0.09 
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